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Abstract. We review the consequences of intrinsic, nonstatistical temperature fluctuations as seen in ob- 
servables measured in high energy collisions. We do this from the point of view of nonextensive statistics 
and Tsallis distributions. Particular attention is paid to multiplicity fluctuations as a first consequence of 
temperature fluctuations, to the equivalence of temperature and volume fluctuations, to the generalized 
thermodynamic fluctuations relations allowing us to compare fluctuations observed in different parts of 
phase space, and to the problem of the relation between Tsallis entropy and Tsallis distributions. We also 
discuss the possible influence of conservation laws on these distributions and provide some examples of 
how one can get them without considering temperature fluctuations. 



PACS. 89.75.-k complex systems - 24.60.-k statistical theory and fluctuations 
duction (relativistic collisions) - 25.75.-q relativistic heavy ions collision 



25.75.Dw particle pro- 



1 Introduction 

Nowadays the statistical approach is a standard procedure 
used to model high energy multiparticle production pro- 
cesses [T]. However, it has been realized that data on many 
single particle distributions deviate in a visibly way from 
what one expects from the usual statistical models, based 
on Boltzman-Gibbs (BG) statistics. These frequently show 
power-like rather than exponential behavior, and, in addi- 
tion, multiparticle distributions are broader than naively 
expected. These observations prompted the idea of a suit- 
able modification of a simple statistical approach used by 
including in it the possibility of accounting for possible 
intrinsic, nonstatistical fluctuations. These were identified 
as the source of the deviations. Such fluctuations are im- 
portant as possible signals of phase transition(s) taking 
place in an hadronizing system [2]. Therefore it is impor- 
tant to be able to include them. In this way the Tsallis 
statistical approach [5], already known in other branches 
of physics, was successfully introduced to the field of mul- 
tiparticle production processefQ- In this approach a new 
parameter, the nonextensivity parameter q appears, which 
is identified with fiuctuations of the parameter T identified 
with the "temperature" of the hadronizing fireball [S]. 

It was shown there that such a situation can only oc- 
cur when the heat bath is not homogeneous and must 



be described by a local temperature, T, fluctuating from 
point to point around some equilibrium value, Tq. Assum- 
ing some simple diffusion picture as being responsible for 
equalization of this temperature [5l|4] one obtains the evo- 
lution of T in the form of a Langevin stochastic equation 
with the distribution of l/T, g{\/T), emerging as a solu- 
tion of the corresponding Fokker-Planck equation. It turns 
out that in this case g{l/T) takes the form of a gamma 
distribution, 
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Convoluting the usual Boltzman-Gibbes exponential fac- 
tor exp(— _E/T) with this g{l/T), one immediately gets a 
Tsallis distribution, hq{E), with a new parameter q, which 
for g — >■ 1, becomes the usual BG distributior(3: 

K{E) = ^^exp^ ( -- ) = 
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^ For details see our previous review [4]. Here we present 
recent developments in this field not covered there. 



Notice that all distributions used here are defined as 
probability density functions with standard normalization, 
J dEhq{E) — 1. This results in the presence of the prefactor 

(2-g)/T. 
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with 



q = i 



where 



Var{T) 



(3) 



(4) 



directly connected to the variance of T. This idea was fur- 
ther developed in 6^ and [7j (where problems connected 
with the notion of temperature in such cases were ad- 
dressed). This forms a basis for so-called superstatistics 
[5]. In what follows, we shall use this approach when dis- 
cussing Tsallis distributions (except of Section H] in which 
we compare it with distribution obtained from Tsallis en- 
tropy). 

It must be mentioned that temperature fluctuations 
(visualized by g > 1 values of the nonextensivity param- 
eter) also allow for a description of the possible energy 
transfer from or to the heat bath [4]. Namely, if Ty is 
a new parameter characterizing such an energy transfer, 
then 

T^T,ff^To + iq-l)n, (5) 

Fig. [T] shows that such an effect is indeed observed [TS]. It 
is caused mainly by the possible energy transfer between 
the central fireball (participants) and nuclear fragments 
passing by without interaction (spectators jj. Notice that 
this energy transfer is only possible in the presence of fluc- 
tuations, i.e., for q > I, when there are no fluctuations and 
q = I one has Tg// = Tq. 



traditional presentation of thermodynamics [IT. In such 
a general approach, the Tsallis distribution ([2]) belongs 
to the class of general admissible distributions which sat- 
isfy thermodynamic consistency conditions and present a 
natural extension of the usual BG canonical distribution 
(I3|). This, together with a recent generalization of classical 
thermodynamics to a nonextensive case presented in |13) . 
form a constructive answer to the critical remarks we en- 
countered concerning the consistency of Tsallis statistics 
with the usual thermodynamics in 1141. 

Applications of Tsallis distributions to multiparticle 
production processes are now numerous. To those quoted 
previously in [4] one should add some new results from 
[TSlfTBlITT] and presented in ^. The most recent applica- 
tions of this approach come from the STAR and PHENIX 
Collaborations at RHIC ^EUdM and from CMS ^, AL- 
ICE [H] and ATLAS P2] Collaborations at LHC (see also 
a recent compilation [53] jf|- In Section [5] we report on new 
results concerning the consequences of temperature fluc- 
tuations in observables measured in high energy collisions 
obtained since our previous review [1] . In Section [3] the 
influence of conservation laws, forcing the use of condi- 
tional probabilities and resulting in g < 1, is discussed. 
In Section m the differences between Tsallis distributions 
as obtained from Tsallis entropy and the concept of su- 
perstatistics is discussed. A possible experimental check is 
proposed. Section[5]is devoted to yet another, not based on 
statistical models, derivation of Tsallis distribution. Sec- 
tion [B] is our summary. 
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Fig. 1. (Color online) Dependence of T^ff on q for different 
energies. RHIC data points are from [ID] whereas NA49 points 
are from [11] (for, respectively, -/s = 6.3, 7.6, 8.8, 12.3, 17.3 
GeV (negative pions). Fits are: Te// = 0.17 - 7.3(g - 1)^ (full 
line, and T^ff — 0.22 — 1.25(g — 1) (dashed line). In both cases 
T^ff is in GeV. 



It is worth mentioning at this point that fluctuation 
phenomena as discussed here can be incorporated into a 



Similar effect is also expected in propagation of cosmic rays 
through the outer space, cf., [9]; we shall not discuss this issue 
here. 



2 Imprints of superstatistic in multiparticle 
processes 

2.1 Multiplicity distributions 

In [29] (cf. also @]) we saw that T fluctuations in the 
form of Eq. ([T]) not only result in power-like behavior of 
single particle distributions, but also in a specific broaden- 
ing of the corresponding multiplicity distributions, P{N), 
which evolve from the poissonian form characteristic of 
BG distributions to the negative binomial (NB) form for 
Tsallis distributions. In short: whenever we have N inde- 
pendently produced secondaries with energies {£"1=1,. ...at} 
taken from the exponential distribution f{E), cf. Eq. ([3]), 
in which case the corresponding joint distribution is given 

by 



fi{Er. 



and whenever 
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* In addition to applications presented in this review, the 
nonextensive approach has also been applied to hydrodynam- 
ical models [24] and to investigations of dense nuclear matter 
[25]. 
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then the corresponding muhiphcity distribution is poisso- 
nian, 
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But whenever in a given process N particles with ener- 
gies {i?i=i,...,Ar} are distributed according to the joint A''- 
particle Tsalhs distribution, 
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(for which the corresponding one particle Tsallis distribu- 
tion function in Eq. ^ is the marginal distribution) , then, 
under the same condition ([7]), the corresponding multiplic- 
ity distribution is the NB distribution [2 7) . 



P(iV) = 



r{N + k) 
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where k ~ 



(10) 



. -1 

For g — >■ 1 one has fe — >■ oo and PU)) becomes a poissonian 
distribution ([5]), whereas for g — >■ 2 one has fc — >■ 1 and 
(fTO|) becomes a geometrical distribution. For large values 
of N and {N) Eq. (|10l) can be written in the following 
scaling form, 



{N)P{N) -ijjiz 



N 



r{k) 



z'=^iexp(-/cz), (11) 



known as Koba-Nielsen-Olesen (KNO) scahng [5{?ll^Bfl . 

Note that, if in the Poisson distribution ([5]) one fluctu- 
ates the mean value, N = E/T (valid for one-dimensional, 
D = 1, case), using its distribution in the form 



9iN) 
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(where g{l/T) is given by Eq. ([T])) then the resulting mul- 
tiplicity distribution 



N 



P{N) = / dNg (N) — exp {-N) 



m 



(13) 



is the NB distribution given by Eq. dTUlrl . 

^ The connection between q and k was first discovered when 
fitting pp data for different energies by means of the Tsallis for- 
mula ([2]) 31 ). The resulting energy dependence of parameter 
q turned out to coincide with that of 1/fe from the respective 
NB distribution fits to the corresponding P{N). It was then 
realized that fluctuations of N in the poissonian distribution 
(O taken in the form of ^p{N / < N >), Eq. HT]), lead to the 
NB distribution (fTD)) . 

® Actually this has been also noted in [32lFf| and recently 
discussed in [33j where the credit in what concerns the origin 
of discussion of such connection between the Poisson and NB 
distributions has been given to [34) . 



2.2 Equivalence of temperature fluctuations and 
volume fluctuations 

The KNO scaling form (fTTj) . with assumed identification 

(where V is the volume of the interaction region) has been 
used in [35 as a starting point for a description of particle 
spectra by means of fluctuations of volume. In this way it 
was hoped to avoid the notion of fluctuating temperature 
discussed here. The results were encouraging. However, for 
constant total energy as assumed in [35^ , E = const, both 
the volume V and temperature T are related, 



this means that 
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and the mean multiplicity in the microcanonical ensemble 
(MCE), N, can be written as 



N = (N) 
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(V) \{T) 



T 
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This implies that both approaches are equivalent and that 
fluctuations of V assumed in |35] in the form given by 
Eq. ([TB| arise as an effect of fluctuations of T considered 
here with g{l/T) given by Ec^. ([1]). This is not assumed 
but derived from the properties of the underlying physical 
process in the nonhomogeneous heat bath. One should also 
remember that UA5 data [3S] show that KNO scaling is 
broken due to the energy dependence of the parameter 
aQ. In fact, as shown in [^, fc"! = -0.104 -f 0.058 In ^i . 
Therefore, in the scenario with fluctuations of the volume 
V, the scaling KNO form of the P{N) used to model these 
fluctuations is a somewhat rough simplification. On the 
contrary, in the scenario of temperature T fiuctuations, 
P{N) is given by a NB distribution, which adequately 
describes the data. 



2.3 Relation between fluctuations observed in different 
parts of phase space 

2.3.1 g-sum rules 

So far, fluctuations of T as introduced in [5] and measured 
by the corresponding parameter q were discussed using 
examples of distributions of longitudinal phase space (in 
the rapidity variable y and integrated over transverse mo- 
menta), dN/dy, and in transverse phase space, dN/dpT- 



^ A possible solution to solve the breakdown of the KNO 
scaling in multiplicity distributions measured in e'^e~ and pp 
collisions has been proposed in [37) . 
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It was found that the corresponding parameters q, q ^ qx 
and q = qL, respectively, are different. Whereas ql — 1 ^ 
0.1 — 0.3 and grows with the energy of colhsion (measured 
mainly in pp and pp collisions), transverse fluctuations are 
much weaker, gy — 1^0. 01 — 0.1 and vary slowly with 
energy (depending only slightly on whether one observes 
elementary collisions or collisions between nuclei) '38i39 . 
As shown in [2^[I] the same fluctuations of T result in 
broadening of multiplicity distributions resulting in its NB 
form as given by Eg. ljlOp . This time the corresponding q 
describes fluctuations in the whole of phase space, with 

^21 



In (38j it was proposed that, because q 



Pt- 



-1 



\T)/{TY 



(i.e., is given by fluctuations of the total temperature T), 
and assuming that CT^(r) = a'^{TL) + a'^{TT), the resulting 
values of q should not be too different from 



(T)2 (r)2 



1. 



(18) 



Therefore, because of the dominance of longitudinal (par- 
tition) temperature over transverse, T^ S> Ty, one should 
expect that q ^ qh- This is indeed observed |38j. This is 
the first sum rule for parameters q obtained from different 
measurements. 

Fluctuations of temperature are usually deduced either 
from data averaged over all other possible fluctuations or 
from data also accounting for fluctuations of other mea- 
sured variables. In this case one can refine the experimen- 
tally evaluated q and, for example, when extracting q from 
distributions of dN/dy, one finds that (cf., [JO] for details) 



def Var{T) Var{z) Var (rriT) 



{TY 



{zY 



[mT 



(19) 



where z = ttit/T (with itit — \/vn? -|-p^). This is the 
second sum rule for the nonextensivity parameters q ob- 
tained from different measurements. It connects total q, 
which can be obtained from an analysis of the NB form 
of the measured multiplicity distributions, P{N), with 
qL — 1 — V ar{z) / {zY , obtained from fitting rapidity dis- 
tributions and Var{raT) /{ititY obtained from data on 
transverse mass distributions. When extracting q from dis- 
tributions of dN/drriT, we proceed analogously but now 
with z — coshy/T. 



2.3.2 Generalized thermodynamic fluctuation relations 

So far, we concentrated only on fluctuations of T. We shall 
continue the discussion by allowing the energy (U), tem- 
perature (T) and multiplicity (N) of the system to fluctu- 
ate and propose to express these fluctuations by the corre- 
sponding parameter q [IT] . Our discussion is based on the 
notion of thermodynamic uncertainty relations discussed 
in [15] . It was suggested there that the temperature T and 
energy U could be regarded as complementary, similarly as 
are energy and time in quantum mechanics. One expects 
from simple dimensional analysis that {k is Boltzman's 
constant) 



Deflnite U (isolation) and definite T (contact with a heat 
bath) to represent the two extreme cases of this comple- 
mentarity. This leads to the so called Lindhard's unce_r_- 
tainty relation between the fiuctuations of U and T 



ujIj + Uj'^ = 



1 



where 



ojl = Var{x)/{xY, (21) 



and this, as was shown in [51], can be generalized to in- 
clude all variables: U, T and A^ by using the nonextensive 
approach. One can then study an ensemble in which the 
energy ([/), temperature (T) and multiplicity (N), can all 
fiuctuate. These fluctuations are then connected by the 
following relation: 



'^N 



1 

Jn) 



- ujrp — 2puJijUJT 



(22) 



{uju - ujt) + 2a;;7WT(l - p) = |(? - 1|, 



where p — p{U,T) £ [—1,1] is the correlation coefficient 
between U and T. This generalizes Linhard's thermody- 
namic uncertainty relation, Eq. ((2T|) . The correlation co- 
efficient enters since when all variables, U, N and T fiuc- 
tuate, the pairs of variables, {U,N) and {U,T), cannot all 
be independent because 

Var{U) = {T)Cov{U, N) + {N)Cov{U, T) (23) 

(cf., [in])- This means that, in general, 

LOU = p(f/, N)u:n + p{U, T)ujT. (24) 

where p{X, Y) denotes the corresponding correlation co- 
efficients between variables X and Y. It should be noticed 
at this point that in the literature [H] there is a similar 
relation connecting the volume, V, pressure, P and tem- 
perature, T: 

UJp = UJy + w^, (25) 

but we shall not discuss it here. 

The observed systematics in energy dependence of the 
parameter q, deduced from presently available data, is 
shown in Fig. [2] From measurements of different observ- 
ables one observes that, for high enough energies, q > 1 
and that values of q found from different observables are 
different. The latter is caused either by technical (method- 
ical) problems or else by a physical cause. The former 
arises when, for example, fiuctuations of temperature are 
deduced either from data averaged over other fluctuations, 
or from more refined data also accounting for fluctuations 
of other variables (as in jlD], see Eq. (|T5|) '). The latter case 
is connected with the fact that the observed g's were ob- 
tained in different parts of phase space. In this case one 
gets an uncertainty relation ((22]) with the help of which 
one can connect fluctuations observed in different parts of 
phase space. For example, one can recalculate q obtained 
from P{N) (i.e., obtained from the whole phase space, see 



AU Al3 > fc, where /3 = 1/T 



(20) 



This idea is still disputable, see jj^j nevertheless we shall 
treat these increments as a measure of fluctuations of the cor- 
responding physical quantities. 
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dashed line in Fig. [5]) and compare it with q evaluated 
from fipr) (i-C, obtained from only transverse part of 
phase space, see full line in Fig. [2]E|. 
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Fig. 2. (Color online) Energy dependencies of the parameters 
q, qL and qr as obtained from different observables. Triangles: 
qL obtained from an analysis of rapidity distributions 31 j solid 
triangles show the uncorrected values, whereas open triangle 
indicates the corrected value [ID]. Squares: q obtained from 
multiplicity distributions P{N) (fitted by g = 1-1- 1/fc with 
1/k = -0.104 -I- 0.029 ln(s)) [17|. Circles: qr obtained from a 
different analysis of transverse momenta distributions, / (pr). 
Data points in this case come, respectively, from the _27, com- 
pilation of data (full symbols) and from CMS data (half filled 
circles at high energies) [20]. The full and dotted lines come 
from Eq. H27p and show, respectively, the energy dependence 
of qr and energy dependence of qL (for p = 0, a — 2/3 and 
K= 1). 
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Fig. 3. (Color online) Example of p obtained from Eq. (|29|) . 
The shaded area shows the extent of possible error, due to the 
uncertainty in fixing gi,- 
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The correlation parameter p appearing here bears im- 
portant information on the details of the production pro- 
cess. For example, p < means that a large energy V (i.e., 
large inelasticity of reaction, K) results in a large num- 
ber of secondaries of lower energies, whereas p > means 
the opposite, one gets a smaller number of larger energies. 
From Eq. (|22[) one finds that the coefficient p is a function 
of all the nonextensivity parameters involved. Denoting by 
Of the part of fluctuations of T in the transverse direction, 
one finds 



qT 



1 — auj, 



T' 



qL-l=ujjj + {I- a)u!' 



and further 



g - 1 = (9L - 1) + (gr - 1) - 2puju^T- 



(26) 



(27) 



^ See |41| for details. A comment is in order concerning re- 
sults of Fig.[2]obtained from f{y). Namely, it turns out that, in 
the fitting procedure, parameters T and q are strongly corre- 
lated [23II40| . As a result q values evaluated in different analysis 
of rapidity distributions [3H] differ slightly from those presented 
here (they give q values comparable or somewhat higher that 
one obtained from multiplicity distribution). 



Fig. 4. (Color online) Dependence of the correlation coefficient 
p on the parameter a for different values of qL ■ 



It can be shown that 

LOT 



V9T- 1 



(28) 



Finally, one obtains correlation coefficient p expressed in 
terms of different fluctuations (in principle measured) (cf. 

ESI " 



m. 



1- 



P = 



(q-i)-(qL-i) 



qr - I 



^V»( 



91.-1 
9r-l 



(jj^ 



(29) 



An example of the feasibility of deducing p from data is 
presented in Fig. [3] for data on p+p at 546 GeV 36. ■ In this 
case one takes from P{N) q — 1.27, from the distribution 



^° In [29) we used a — 2/3 and k = 1; for p — 0. However, the 
actual values of a and k parameters are irrelevant in this case. 
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oi pt one has qx = 1.09, whereas from the original ql = 
1.36 one obtains, after correction, q^ — 1.14 (cf. Fig. [2]). 

To summarize this part, note that, to get the correla- 
tion coefficient /9, one has to know all the fluctuations, i.e., 
both in the entire phase space, q, as separately in its trans- 
verse, qr, and longitudinal, qL, parts. The best known is 
q (no corrections needed) , for qx the corrections are small 
and can be neglected, finally, for q^ the corrections are 
large and must be accounted for (cf.. Fig. [3]). 



2.4 Energy fluctuations - heat capacity 

We now present energy fluctuations resulting from Tsal- 
lis statistics and emerging from our analysis [151[^ . This 
subject already has its history (cf. [47]) and was also re- 
cently under investigation (cf. [i5]). 

In Boltzman statistics [551 (with kT = 1/13 — const 

and A'^ = const) the energy U — J2i=i ^i ^^ ^ particles 
is distributed according to 



9t,n 



/3 



r{N) 



m)"''^M-m 



(30) 



for which 
VarjU) 



MB) N 



^ where Cy = ^W-- (31) 



dT 



In Tsallis statistics [15] one has, respectively, 
■ [l-(l-9)/3[/]T^+'"^ 



for which 



Var{U) 
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TV 



39 \C^P 
9-1 
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(32) 



(33) 



39 
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where 



C. 
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d{U) 
dT 



Nk- 



3-2q 



a 



(B) 



3-2q 



(34) 



Notice that fluctuations of the energy U are, in general, 
given by the sum of two components: one obtained in the 
case of no fluctuations and given by the heat capacity Cy ' 
(which we call the kinetic component) and one originating 
in fluctuations, and given by the heat capacity C^f' (van- 
ishing when fluctuations vanish, we call it the potential 
component): 



k 



k 



VarjU) 



(35) 



where (cf. [47 ) 

4^)=yfcTV and C^f) ^ k^^^^^. (36) 

^ N +\ q-\ ^ ' 

From analysis of nuclear collisions we know [15_ that q 
depends on N: 



9-1 = 



a 

iV' 



(37) 



where a is some constant of order unity depending on the 
reaction considered. We can therefore write 
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4-3g 
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(38) 



For small values of q — 1 (in practice already for q— 1 << 
0.5) one has 



Var{Ul 1 

^^7y^>('^-i)- 



1 + a 

N ■ 



(39) 



3 Conditional probability - influence of 
conservation laws 

Let {-El,..., at} be a set of N independent identically dis- 
tributed random variables described by some parameter A 
and let g]^{E, A) denote the gamma density function with 
parameters N and A. For independent energies, {£"2=1,..., at}, 
each distributed according to the simple Boltzman distri- 
bution: 

1 / E,: 



51 iE^) 



A 



the sum 



exp 
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A 



E = Y^E, 



(40) 



(41) 



i=l 



is then distributed according to the following gamma dis- 
tribution, 



9n{E) 
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A(TV-l)! VA 
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E 
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(42) 



If the available energy is limited, for example ii E = 
J2i=i Ei = Na = const, then we have the following con- 
ditional probability for the single particle distribution, 

f{E^)■. 



f{E,\E^Na) 



gi{E,)gN-i{Na-E,) 



gN{Na) 



jN-l) ! 
N a 



1- 



N a 



N~2 



(43) 



This is nothing else then the well known Tsallis distribu- 
tion 



f {E,\E ^ const) = 



2-q' 



l-(l-9')f 



(44) 
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with 



N ~3 



q = 



< 1 and A = (3 - 2q') a 



(45) 



N -2 

which is always less than unity. Here A — const and do 
not fluctuate. 

Now consider a situation in which the parameter A in 
the joint probability distribution 

N 

9{{El.^...,N})^l[9r{Ei) 

i=\ 

fluctuates according to a Gamma distribution, Eq. ([T]). In 
this case we have the single particle Tsallis distribution 



which, when additionally iV >> 1 (or g' ^- 1) reduces to 
Eq. ge]). 



h^ (E,) 



1-(1"?) 



E, 



(46) 



10' 
10-^ 
10-' 
10 = 

S 10"' 

Q. 

10-^ 

10" 
10" 



and the distribution oi E = J2i=i Ei is given by (cf. [22]): 



10" 



N 


w^j^^ 












^. 






\ 


•\1 




q'=0.964 (N=30) 


\ 




■i 


q=1.0 




• ■ 


n ■ 
' 1 


q=1.1 
q=1.2 




, 1 I 1 ' 


• 

• ■ 

1 1 1 : T 



hNiE) 



Ar(7V)r(^) 



7V-1 



i-(i-,)f 



1-^+T^ 



(47) 



0.0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1.0 
X=E/E 



Fig. 5. (Color online) Conditional probability distribution, 
P{x = E^/E), for g = 1 ( Eq.gH)) and g > 1 (Eq.^, in 
both cases TV = 30 (q' — 0.964), compared to exponential dis- 
tribution {q = 1, q' — 1).). 



^N 



If the energy is limited, i.e., ii E = J2i=i Ei — Na = 
sonst, we 

f{E^\E) 



const, we have the following conditional probability: 

h,{E,)hN-i{E-E,) 



hNiE) 
{N^l){2-q) 



A' fE-E, 



E[{3-2q)-N{l-q)] X \ E 



N-l 



l-(l-9) 



l + (l-,)f 



1-g 

(48) 



where 

A' = A - (1 - q)E. (49) 

For gr ^ 1 Eq. (gH]) reduces to Eq. ^. On the other 
hand, for large energy (E — >■ oo) and large multiplicity 
{N — >■ oo), the conditional probability distribution (|^5)) re- 
duces to the single particle distribution given by Ea. (|46)) . 
Introducing the parameter q' defined in Ea.(|45p the con- 
ditional probability ([^H)) can be rewritten as 

(2 - q') (2 - q) A' 




f{Em 



E[{3 - 2q) {1 - q') - {3 - 2q') {1 - q)] A 
E-E,\^ 



E 



!-(!-, 4 



E, 



i + a-'z)Tf 



A' 



(50) 



For Ei << E it becomes 

{2 - q') {2q - l){q - 1) 



0.0001 



Fig. 6. (Color online) Ratio of conditional distribution func- 
tion f{Ei\E) and single particle distribution hi{Ei) as function 
of x = Ei/E for Tsallis statistics {q = 1.1 and N ^30). 



The results presented here are summarized in Figs. [S] 
and|6]which shows how large differences are (in x — Ei/E) 
between the conditional Tsallis distribution f{Ei\E) and 
the usual fei ('£',■ r^. 



fiEr\E) 



A[(3 - 2q) (1 - q') - (3 - 2q') (1 - q)] 

1 

1-g 






(51) 



^^ We would like to stress that Ea. H43|l has the form of a mi- 
crocanonical distribution in the one dimensional case, D = 1. 
In [49] it was shown that smearing this distribution over a 
Gamma type multiplicity distribution results in a microcanon- 
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4 Tsallis entropy and the Tsallis distribution 
function - nonadditivity in nuclear collisions 

In all examples discussed so far we treated the Tsallis dis- 
tribution, Eq. ([5]), as a kind of superstatistics [5] without 
really resorting to Tsallis entropy ^. However, closer in- 
spection of both approaches reveals that the correspond- 
ing nonextensivity parameters (say q and q' , respectively) 
are not identical. In fact one encounters a sort of dual- 
ity, like q = 2 — q' discussed, for example, in [STII^FT] . We 
shall now address this problem in more detail (cf., [T7] for 
details) . 

When starting from Tsallis entropy |3J, 



Sa 



1 



1-g 



dxf'ix) - 1 



(52) 



one can obtain the probability density function f{x) cither 
by optimizing it with constraints 



(53) 



dxf{x) = 1; / dxxP{x) = {x)g, 

in which case [55] 

fix) = {2-q)[l-{l-q)x]^; 

0<a;<oo; l<q< 3/2, 

or else by using as constraints 

dxf{x) = 1; / dxxf{x) — (x). 



(54) 



in which case 1521 



(55) 



(56) 



[1 + (1 - q)x] 1-' ' 

< a; < oo; 1/2 < q < 1. (57) 



/(^•) = 



Out of these two possibilities, only (IMI) is the same as the 
distribution obtained in superstatistics and used above, 
cf., Eq. (I2|). On the other hand, the second distribution, 
Eq. (j56p . which seems to be more natural from the point 
of view of a physical interpretation of the constraint used, 
becomes the first one if expressed in terms of q' given by 



Namely, in this case one has 



fix) = i2-q')[l-il-q')x]- 



(58) 



(59) 



which, as show in Fig.[2j when compared to single particle 
distributions, results in q' > 1. 



leal generalization of the Tsallis distribution which fits the 
fragmentation functions measured in e'^e~ experiments with 
similar q{s) evolution to that presented in Fig.[51 It was demon- 
strated that this type of energy dependence seems to be con- 
sistent with the DGLAP evolution equations |50] . 



It turns out that there are data allowing the above du- 
ality (at least in principle, considering the present status 
of their quality). They are provided by nuclear collisions 
in which one observes the apparent nonadditivities which, 
as will be shown, allow us to compare and discuss both q 
and gH 

We start with the phenomenological approach used to 
describe nuclear collisions which is based on the superposi- 
tion model with main ingredients being nucleons that have 
interacted at least once [53] ■ In this case, when sources are 
identical and independent of each other, the total (N) and 
the mean {{N)) multiplicities are supposed to be given by. 



iV.^. 



and 



{N) = {,y){n, 



(60) 



where v denotes the number of sources and rii the multi- 
plicity of secondaries from the i*'' source. Albeit at present 
nuclear collisions are mostly described by different kinds 
of statistical models [1 , which automatically account for 
possible collective effects, nevertheless a surprisingly large 
amount of data can still be described by assuming the 
above superposition of independent nucleon-nucleon colli- 
sions (possibly slightly modified) as the main mechanism 
for the production of secondaries. The question of the 
range of its validity is a legitimate one [Slj . 

Using the notion of entropy, and considering v inde- 
pendent systems for which the corresponding individual 
probabilities are combined as 



P^^^ {^1 



Wp'iH^k), 



(61) 



k=l 



and assuming that all pq (xk) are the same for all k (i.e., 
their corresponding entropies 5*^ are equal), one finds 



4^^-E7^a-^r1^f^ 






-1 



Notice that 




In 


1 + (1 - 


- q)Si'''> 


and that 








S]^ 



— v\n 



q^l 



1 + (1 - g)5W 



s\ 



(1) 



For q < 1 one has 



S, 



(i') 



oo, 



v 



(62) 



(63) 



(64) 



(65) 



^^ Apparently similar duality occurs in nonextensive treat- 
ment of fermions for which the particle-hole correspondence, 
nq{E,T,fi) = 1 — n2-q{—E,T,—ii) (where fi is the chemical 
potential), must be preserved by the q- Fermi distributions |25) . 
However, here we are facing different problem, namely that pa- 
rameter q in entropy Sq differs from parameter q' in probability 
distribution fgi with q — 2 — q' . 
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i.e., entropy Sq is nonextensive. For q > 1 one has 
S"^") > only for g < 1 + -\- 



On 



(66) 



and 



S, 



H 



0, 



i.e., entropy is extensive. 



0<f^<5(i). 



(68) 



In the foUowing we put v = Nw/2 = Np {Nw is the 
number of wounded nucleons and Np is the number of 
participants from a projectile). Assuming naively that the 
total entropy is proportional to the mean number of pro- 
duced particles, 

S = a{N), (69) 

one obtains the following relation between mean multi- 
plicities in AA and NN collisions, 



a{N) 



AA 



[l + {l-q)a{N) 



^Np 



PPi 



1 



l~q 



(70) 



At this point we stress the following observation, so far not 
discussed in detail. Namely, because (as shown in [T5]). 
{N)aa increases nonlinear ly with Np and {N)aa > Np ■ 
{N)pp, the nonextensivity parameter obtained here from 
considering the corresponding entropies must be smaller 
than unity, q < 1. On the other hand, all estimations of the 
nonextensivity parameter (let us denote it by q') discussed 
before lead to g' > 1. This is the q duality in nonextensive 
statistics mentioned above, on which we shall concentrate 
in more detail. 

To start with, the relation (Tzni) is not exactly correct 
for Sq. In what follows we denote entropy on the level of 
particle production by s (and the corresponding nonex- 
tensivity parameter by g), whereas the corresponding en- 
tropies and nonextensivity parameter on the level of A^ A^ 
collisions by S and q. From Eq. ([S5)) we have that for A^ 
particles 



*9 



1 + (1-9)4 



(1) 



N 



1 



1-? 



"^N-s^}^^aN, (71) 



„(i) 



where s~ = a is the entropy of a single particle. In a 
A + A collision with v nucleons participating, Eq. (|65l) 
results in 



'=>q - 



l + {l-q)S'q 



(1) 



- 1 



1 



(72) 



:(i) 



where Sq is the entropy of a single nucleon. 

Denoting multiplicity in a single N + N collision by n, 
the respective entropy is 



^q — ^q 



(1) 



1 + (1 - q) s 



(1) 



-1 



(73) 



whereas the entropy in a A + A collision for N produced 
particles is 



(67) This 



q{N) 

^q 



that 



l + {l-q)s, 



(1) 



N 



-1 



l-q 



"^q — '^q ■ 



(74) 



(75) 



Notice that parameters q and q are usually not identical. 
Moreover, from the relation 



q-l = 



1 



aNp 



, Np\ Cv 



A 



Np 



(76) 



one finds that for A^A^ collisions (where Np — A) q — 1. 
On the other hand, for g = g Eq. ((75)) corresponds to the 
situation encountered in superpositions, as in this case one 
has 



1 + (1-9)4'^ 



N 



1 + (1-9)4^ 



and so 



N = nv. 
Consider now the general case and denote 



(1). 



ci = 1 + (1 - g) 4^^ 



1 



C2 



1 



These quantities are not independent because: 

C2cf + l-C2 = {C2C1 + 1 - C2Y . 

From relation (150)) 



1 



A^ 



V ■ n vn Anci 



In 



{C2C^ + 1 - C2T - (1 - C2) 



C2 



(77) 
(78) 

(79) 
(80) 
(81) 



which for A^ = {Naa), n — (Npp) and v = Np is presented 
in Fig. [7| for different reactions. As seen there one can 
describe experimental data by using C2 = 1.7 and with 
ci depending on energy ^/s according to ci(s) — 1.0006 — 
0.036s^^°'^^. Notice that for energies ^/s > 7 GeV one has 
ci > 1. This means that q < 1 and (because C2 > 0) also 

q<l- 

To summarize this section, we have shown that, non 
additivity in the superposition model described using the 
notion of entropy clearly requires q < 1, cf. Figs. [7| This 
means that q' is not the same as q. The conclusion one can 
derive from these considerations is that the second way of 
deriving /(a;), which uses a linear condition, cf. Eq. ()56|) . 
is the correct one, and that q' in the distribution is not 
the same as q in the entropy. The problem is that, whereas 
from distributions one can easily deduce a numerical value 
of q' , this is not the case when one uses entropy (at least 
not when deduced from presently available data). There 
are too many variables to play with (cf., considerations 
using the superposition model as above). For example, in 

the definition of ci in Eq. (|79)) . one has s- , which is not 
known a priori. The only thing one can deduce in this 
case is that q < 1. We cannot therefore check numerically 
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Fig. 7. (Color online) Energy dependence of the charged mul- 
tiplicity for nucleus-nucleus collisions divided by the superpo- 
sition of multiplicities from proton-proton collisions (cf. Eq. 
((ST])). Experimental data on multiplicity are taken from the 
compilation |55] . 



on Tsallis entropy, are by no means the only possibili- 
ties. Therefore we end with short discussions of two ex- 
amples of obtaining Eq. ^ in a completely nonthermal 
way, these are the application of order statistics and the 
use of stochastic networks. 



5.1 Order statistics 

Order statistics is based on the observation |S3] that the 
selection of the minimal value of the ordered variables 
leads in a natural way to its distribution being given Eq. 
(I2|) (with q both greater and smaller than unity, depend- 
ing on circumstances), i.e., in fact by the Tsallis distribu- 
tion , the same as that resulting from Tsallis nonextensive 
statistics. Distribution of the minimal values of some spe- 
cific choices of the variable E is known in the literature as 
order statistics 64 O- 



that relation ([58)) really holds. But, if one agrees that the 
Tsallis distribution comes from Tsallis entropy, we have 
only two options: either q' = q or q' — 1 = I — q. Our 
conclusion presented here, that q' > 1 and q < 1, therefore 
supports the second option, i.e., Eq. ((58|) . 

However, this final observation calls for comment. Namely, 
the probability density function (PDF) is usually evalu- 
ated by the Maximum Entropy Method (MEM) for Tsal- 
lis entropy with some constraints [S5^ Ej- Therefore the 
situation is not unique since there are four possible well 
documented MEMs [SJ using two kinds of definition for 
an expectation value of the physical quantities: the nor- 
mal average (1551) and the ^-average (|53)) (with normal, as 
here, or the so-called escort PDFs ^H])- Although various 
arguments justifying it have been given [59i it was also 
been pointed out that, for a small change of the PDF, 
thermodynamic averages obtained by the ^-averages are 
unstable, whereas those obtained by the normal average 
are stable [60] . On yet another hand, it is claimed that for 
the escort PDF, the Tsallis entropy and thermodynamical 
averages are robust [ST] . All this means that the stability 
(robustness) of thermodynamical averages as well as the 
Tsallis entropy is still a controversial issue ^5] . 




100 



1000 



Fig. 8. (Color online) q as function of n given by Eq. (|84 
different values of n. 



for 



5 Examples of nonfluctuating (nonthermal) 
mechanisms leading to Tsallis distribution 

It should be realized that the so far discussed origins of 
the Tsallis distribution, based either on superstatistics or 



We now present a generalized version of what was pro- 
posed in [53] . We start with a set of n virtual particles (so 
called ghost-particles) with energies Si taken from some 
distribution /(e). Ordering the values of Si (i.e., intro- 
ducing in this set rank statistics), £i < £2 < . . . < £«, 
we choose a real particle with minimal energy E ^ ei = 
min {{ei}). It is straightforward to find a function g{E) de- 
scribing the energy distribution of real particles. The prob- 
ability density to find a particle with energy E among n el- 



Notice that Tsallis entropy is a monotonic function of the 
Renyi entropy, Sq = In^ [exp (Rq)], and both lead to the same 
equilibrium statistics of particles (with coinciding maxima in 
equilibrium for similar constraints on the expectation value). 



^ Actually, one can easily invent a nonthermal scenario lead- 
ing to a thermal-like form of the observed spectra, see, for 
example, recent work [65j . In such an approach the resultant 
distribution emerges not because of the equilibration of ener- 
gies due to some collisions (i.e., because of the kinematic ther- 
malization), but rather because of the process of erasing of 
memory of the initial state and is the result of the approach- 
ing to a state of maximal entropy (called in [65] stochastic 
thermalization) . 
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ements is nf{E). The probability to find particles with en- 
ergy greater than E is 1 — F{E), where F{E) = L def{e) 
is the distribuant of / If a particle of energy E is already 
that of the minimal energy it means that the remaining 
n— 1 particles have to poses higher energies. The probabil- 
ity of such an event is equal to [1 — F{E)]^~-^. This means 
that the distribution of the minimal value in sample of n 
elements ij^ 



giE)^nf{E)[l~-FiE)Y 



(82) 



Because f{E) — dF{E), the distribution g{E) is properly 
normalized if f{E) is normalized. For 



f{e)^-a{K + l){l + aey 



(83) 



where k 7^ — 1 (because of the normalization requirement) 
and a = —sign{K + l)/3 (/3 — l/T > 0) one gets g{E) in 
the form of Tsallis distribution, Eq. ^ , with 



q = 



n(K + l)- 2 
n(K + l)- 1' 



(84) 



q > I ioY K < —1 and q < 1 ior k > —1. Fig. |H]shows q{n) 
dependence for different values of the parameter k (special 
cases of K = — 2 and k = were discussed in [55]). 



5.2 Stochastic networks 

Stochastic network structures occur in almost all branches 
of modern science (including sociology and economy). They 
have therefore been the subject of intensive research, also 
by means of Tsallis statistics (cf. [66ll67] for details and 
full list of references; in fSS" this approach has been ap- 
plied to multiparticle production processeto). There are 
two basic types of stochastic networks: 

— Networks with a constant number of nods, M, for 
which probability that given node has k connections 
with other nodes (fc links) is poissonian 72 , 



i^(fc) = f ---; 



kq = (k). 



(85) 



Networks in which the number of nodes is not station- 
ary and the distribution of links P{k) is given by dy- 
namics of the growth of network [73^ . It varies between 
being exponential, 



P{k) ~ — ■ exp 



(86) 



^^ More formally, the cumulative distribution function is 
G{E) = 1 - [1 - F{E)]" and the density distribution is 
g{E) = dG{E)/dE = n[l - F{E)]"-^dF{E)/dE = nf{E)[l - 
F{E)]"-'^. 

^^ In [68] the "power laws", assumed ad hoc in [69] (as a 
kind of opposition to Tsallis statistics), was explained using 
a stochastic networks approach presented here. Actually, this 
"power laws" idea is continued recently in |70) as an appar- 
ent new observation. It must be mentioned therefore that this 
idea is actually quite old; such a type of parametrization of 
Pt distributions has been proposed (and was shown to be phe- 
nomenological successful) already in [71) . 



and power-like, 



P{k) 



2kH 

Kq + t 



(87) 



behavior. In the former case each new node connects 
with the already existing ones with equal probability, 
n{ki) = 1/(ko -\-t — 1), independent of fc^. In the lat- 
ter case one has preferential attachment (the so called 
" rich-get- richer" mechanism, here k, < kq is the num- 
ber of new nodes added in each time step) with, in this 
case, n{ki) = ki/{2Kt) choice. 

Let us remind ourselves that, whereas 



ax A 



/(a;) == ^ cxp 



for the x-dependent scale parameter 

A —> X{x) = Ao - (g - l)x 
the exponential solution takes a power-like form, 

2-q 



(89) 



/(^) 



An 



1 



(l-9)f 
Ao 



For preferential attachment used in ( 
ter equation" 

-^r = -'^^(^^ 

by the assumed " growth of the network" 

-— = a + bk, 
at 



(90) 

dividing the " mas- 

(91) 

(92) 



one obtains the following evolution equation for the net- 
work considered: 



dP{k) 
dk 



-cP{k) 



dt_ 
'dk 



bk 



P{k). 



(93) 



For c = 1, a = Ko ^nd b = q — 1 one has K(fc) = kq + [q — 
l)k and a solution of Eq. ([M]) in the form of the Tsallis 
distribution, Eq. ([^0)) : 



P(fc) 



KO 



l-(l-g)- 



KO 



(94) 



For q — > 1 Eq. ([M]) recovers Eq. (|551) whereas for k » 
KQ/{q — 1) it leads to "scale- free" power distribution 



Pq{k) cc k '^ , with 7 



1 



1 



(95) 



The frequently observed value 7 = 3 therefore corresponds 
to q = 4/3. At this value of q the variance of distribution 
P{k) diverges, 



Var{k) 



nl{2-q) 



9^4/3 



(3-2g)2(4-3g) 



(96) 
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We close this section by noticing that formally we can 
interpret Eq. (l93l) as the stationary solution, of the follow- 
ing Fokker-Planck equation, 



d{K2P{k)) 
dk 



KiPik), 



(97) 



where Ki = q—2 and K2 = Ko + {q—l)k, This corresponds 
(cf. network growth given by Eq. (IMjl ) to the Langevin 
equation with multiplicative noise (ry) in the form |S]: 



dk 
'di 



T]k^£„ 



where ^ is the traditional noise term. In this case both 
noises have nonzero mean values: (^(t)) = Ho and {Ti{t)) = 
1 — g, and correlations: Cov {^{t),^ {t')) — 2koS {t — t'), 
Cov {rj{t), 77 (i')) = and Cov (r;(i), ^ (i')) = {l-q)S {t - t'). 



6 Summary 



The possibility of occurrence of intrinsic, nonstatistical 

temperature fluctuations has far-reaching consequences which Phys. E 16, 1687 (2007) 
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